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By  A.  P.  Basu  and  Nader  Ebrahimi 
University  of  Missouri,  Columbia 

ABSTRACT 

Various  definitions  of  multivariate  harmonic  new  better  than 
used  in  expectation  (MHNBUE)  life  distributions  are  introduced  and 
their  interrelationship  is  studied.  These  are  multivariate 
generalizations  of  the  largest  available  univariate  class  of 
distributions  with  aging  properties.  Examples  are  given  to  illus¬ 
trate  these  concepts.  -Various  closure  properties  of  MHNBUE  distri¬ 
butions  are  proved. 


Key  Words  and  Phrases:  Multivariate  distribution,  Multivariate 

aging.  Reliability,  Distribution  with 
aging  property. 


•This  research  was  supported  by  the  NSF  Grant  TM7  8009463 


and  ONR  Grant  N000 14-78-C-0655 . 


Multivariate  Harmonic  New  Better  Than 
Used  in  Expectation  distributions* 


By  A.  P.  Basu  and  Nader  Ebrahimi 
University  of  Missouri,  Columbia 

1.  Introduction.  In  reliability  theory  various  concepts  of 
(univariate)  aging  or  wearout  have  been  proposed  to  study  lifetimes 
of  systems  and  components.  The  five  most  commonly  studied  class  of 
distributions  are  the  following: 

1)  The  increasing  failure  rate  class  (IFR) ;  2)  the 

increasing  failure  rate  average  class  (IFRA) ;  3)  the  new 

better  than  used  class  (NBU);  4)  decreasing  mean  residual  class 
(DMRL) ;  and  5)  the  new  better  than  used  in  expectation  (NBUE) 
class.  For  a  description  of  some  of  these  classes  see  Barlow  and 
Proschan  (1975) .  Recently  Rolski  (1975)  proposed  a  new  class  of 
distributions  called  the  harmonic  new  better  than  used  in  expecta¬ 
tion  (HNBUE)  class  which  will  be  defined  later.  Each  of  the  above 
six  classes  have  their  dual  with  standard  nomenclature.  The  dual 
of  HNBUE  class  is  said  to  be  harmonic  new  worse  than  used  in  expec¬ 
tation  (HNWUE) .  Klefsjo  (1980)  has  studied  the  properties  of 
HNBUE  (HNWUE)  classes  of  distributions.  He  has  proven  several 
closure  theorems  for  this  class  and  the  following  chain  of  implica¬ 
tion  exists  among  the  six  classes  of  distributions. 


♦This  research  was  supported  by  the  NSF  Grant  INT  8009463 


and  ONR  Grant  N00014-78-C-0655 . 
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(DFRA) 
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Thus  HNBUE  is  the  largest  available  class  of  distributions  with 
aging  property. 

Recently  attention  has  been  directed  towards  extending  the 
concepts  of  univariate  aging  to  the  multivariate  case,  (see  for 
example  Harris  (1970),  Basu  (1971),  and  Brindly  and 
Thompson  (1972)),  and  using  those  multivariate  concepts  to  define 
corresponding  multivariate  class  of  distributions.  Buchanan  and 
Singpurwalla  (1977) ,  Block  and  Savits  (1980) ,  and  Ghosh  and  Ebrahimi 
(1980)  have  considered  the  cases  for  multivariate  IFR,  IFRA,  NBU, 

DMRL,  and  NBUE  distributions.  The  purpose  of  this  note  is  to 
propose  multivariate  versions  of  HNBUE  distributions  (MHNBUE) .  It 
is  shown  that  these  are  the  largest  available  classes,  and  include  the 
class  of  MHNBUE  proposed  by  Klefsjo  (1980)  for  the  bivariate  case. 

In  section  2  of  this  paper,  we  have  introduced  the  various 
definitions  of  the  MHNBUE  involving  a  certain  hierarchy,  and  have 
described  their  physical  implications.  Our  definitions  of  the 
MHNBUE  are  different  from  Klefsjo' s  definitions  of  MHNBUE.  We  have 
compared  our  definitions  of  MHNBUE  with  Klefsjo  definitions  of 
MHNBUE.  We  have  also  examined  in  this  section  how  several  impor¬ 
tant  classes  of  life  distributions  satisfy  one  or  the  other  defini¬ 


tion  of  MHNBUE. 
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The  dual  of  the  HNBUE  is  HNWUE.  In  section  2,  multivariate 
HNWUE  (MHNWUE)  definitions  are  also  given  parallel  to  those  for 
the  MHNBUE. 

Various  closure  properties  of  MHNBUE  distributions  under 
different  definitions  are  studied  in  section  3.  It  is  known  that 
in  the  univariate  case  NBUE  class  is  included  in  the  HNBUE  class. 

It  is  examined  in  section  4  how  far  the  MHNBUE  distributions  as 
introduced  by  Buchanan  and  Singpurwalla  (1977)  and  discussed  in 
more  detail  by  Ghosh  and  Ebrahimi  (1980)  lead  to  one  or  the  other 
MHNBUE  definition  as  introduced  in  section  2. 

2.  MHNBUE:  definitions  and  example.  Let  X^,...,Xp  denote  the 
survival  (failure)  times  of  p  devices  having  a  joint  distribution 
function  H  (x.,...,x  ).  The  joint  survival  function  of  these  p 

r  r 

devices  is  denoted  by  Hp(x^, . . . ,xp)  *  PtX^  >  x^,...,Xp  >  xp) .  It 
is  assumed  that  Hp(0,...,0)  *  1.  In  the  univariate  case,  a 
non-negative  random  variable  X^  is  said  to  have  a  HNBUE  (HNWUE) 
distribution  i f 

f"  - 

Jt  H^(x)dx  <.  (_>)y  exp  (-  t-,/y)  for  all  t^  >_  0,  (2.1) 

where  it  is  assumed  that  [  H.(x)dx  =  y  <  «. 

J0 

With  the  alternate  representation  of  (2.1) 


(2.2) 
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H.(y)dy 

where  e„  (x)  =  ■-» — — - ,  it  is  easy  to  see  that  the  condition 

H1  Hx  (x) 

is  equivalent  to  saying  that  the  integral  harmonic  mean  value  of 
the  mean  residual  life  of  a  unit  at  age  x  is  less  (greater)  than 
or  equal  to  the  integral  harmonic  mean  value  of  a  new  unit. 

Definitions.  Our  definitions  of  MHNBUE  are  the  natural  multivariate 
extension  of  (2.2).  Suppose  for  simplicity  that  H  (x-,...,x  )  >  0 

r  A  r 

and  let 


V*1 


•v  * 


•00  0<X> 

JXj' '  ’  I  x  Vyl . Vdyl"'dyl 


H  (x.. .... i x  ) 

pi  p 


(2.3) 


(x. 


,x) 


I  H  (y,...,y)dy 
V 


Hp  ( x  j  .*. .  ,x) 


(2.4) 


H  is  said  to  be 
P 


(i)  MHNBUE-I  (MHNWUE-I)  if 


P 


eH1(xl' *  *  * 'xp>dxi'  dx2 '  *  *  * 'dxp 

P 


< 


(i)en  (0 


0)  (2.5) 


for  all  t^iO  (1  <  i  <  p),  and  similar  inequalities  are  assumed  to 
hold  for  all  subsets  of  random  variables. 
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(ii)  MHNBUE-II  (MHNWUE-II)  if 

tP 

- <  (>)eH  (0,...,0)  (2.6) 

(x« « . . » r x  ) dx. ... dx  p 

1  pip 

for  all  t  >.  0,  and  similar  inequalities  hold  for  all  subsets  of 
random  variables. 

(iii)  MHNBUE-II I  (MHNWUE-III)  if 

t 

“7t  “ 

e  (x,...,x)dx 

'0  p 

for  all  t  >.  0,  and  similar  inequalities  hold  for  all  subsets  of 
random  variables. 

(iv)  MHNBUE-IV  (MHNWUE-IV)  if 

<.  gH  (0,  .  .  .  ,0)  (2.8) 

P 

for  all  t  >.  0,  and  similar  inequalities  hold  for  all  subsets  of 
random  variables. 

interpretation  of  these  four  defini- 

P 

IT  (y.  -  x.  )dH(y.  , .  . .  ,y  )  ,  (2.9) 

i=l  111  p 

for  proof  see  Ghosh  and  Ebrahimi  (1980).  Now  using  (2.9)  it  follows 


Next  we  give  the  physical 
tions.  First  note  that 


!h  (xi . V  ‘I  -I 

p  1  ;xn  Jy 


f  g„1(x, . . .  ,x)dx 

'0  Hp 


.  (>)eH  (0, 


,0) 


(2.7) 


that  (2.5)  is  equivalent  to  the  statement  that'  the  integral 
harmonic  mean  value  of  the  conditional  mean  residual  product  life¬ 
time  of  the  components  of  a  unit  with  the  components  surviving 
ages  x^,...,Xp  respectively  is  less  (greater)  than  or  equal  to 
the  integral  harmonic  mean  value  of  the  mean  product  lifetime  of 
the  components  of  a  new  unit.  Similar  interpretation  holds  for 
(2.6).  The  definition  (2.7)  is  equivalent  to  the  statement  that 
the  integral  harmonic  mean  value  of  the  conditional  mean  residual 
product  lifetime  of  the  components  of  a  unit  when  all  the  components 
have  survived  a  certain  time  X  is  less  (greater)  than  or  equal  to 
the  harmonic  mean  value  of  the  mean  product  lifetime  of  the  compo¬ 
nents  of  a  new  unit.  Finally,  definition  (2.8)  is  equivalent  to 
the  statement  that  a  multivariate  distribution  is  HNBUE  (HNWUE)  if 
the  minimum  of  the  components  has  a  univariate  HNBUE  distribution. 

It  is  trivial  to  check  that  MHNBUE-I  —  =>  MHNBUE-II.  However, 
the  following  examples  show  that 

MHNBUE-I  =t^=>  MHNBUE-IV  (so  that  MHNBUE-II  =£*  MHNBUE-IV)  . 

Example  1.  Let  and  iid  with  common  survival  function 

II  if  0  <  x  <  3 
1/4  if  3  <  x  <  7  . 

0  if  x  >  7 

Then  by  using  (P2)  of  section  3,  (X^,  X2)  is  MHNBUE-I.  But 
min^,  X2)  is  not  1INBUE»  i«e*  •  (X^  *2)  is  n0t  MHNBUE-IV* 

The  next  example  shows  that  MHNBUE-I  =7*=*  MHNBUE-III  (so  that 
MHNBUE-I  I  =/=*  MHNBUE-IV)  . 
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Example  2.  Let  X1  and  X2  be  independent  with  survival  functions 


Fv  (x) 
A1 


1  if  0  <  x  <  3 
1/4  if  3  <  x  <  7 
0  if  x  >.  7 


and  F  (x) 
X2 


1  if  0  <  x  < 
1/2  if  1  <  x  < 


1 

5. 


0  if  x  >.  5 


respectively.  By  using  (P2)  section  3,  (X^,  X^)  is  MHNBVE-I. 

But  (Xx,  X2)  is  not  MHNBUE-III. 

Finally,  the  following  example  shows  that 
MHNBUE-IV  =/=>  MHNBUE-III. 

Example  3.  Let  X^,  X2  denote  the  survival  times  of  2  devices  having 
a  joint  survival  distribution  function 

H2(x1,  x*2)  =  exp(-  max(X3x1,  X4x2) ) , 

where  A.J,  X4  >  0.  Then,  (X1#  X2)  is  MHNBUE-IV.  But  (X1#  X2)  is  not 
MHNBUE-III. 

Remark  1.  Example  1  shows  that  MHNBUE-III  MHNBUE-IV. 

Klefsjo's  (1980)  definitions  of  tjie  MHNBUE  are  the  natural 
multivariate  extension  of  (2.1). 

We  first  introduce  the  following  definitions  of  MHNBUE  as 
given  in  Kief es jo  (1980) •  A  bivariate  distribution  function  H2  is 
said  to  be  MHNBUE  if 


0OO  f  00  00 

Hj  (Xf  +  t^,  x2  +  t2)dt1dt2  <  j  GfXj^  +  tj^ ,  x2  +  t2)dt1dt2 


for  x^»  x2  >.  0; 
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B. 


r  "2ui 


+  t,  x2  +  t)dt 


-  1  GUi + 1- 


x2  +  t) dt 


for  Xj.,  x2  >  0; 


n00  #00  #00 

_  „  ff2<x  +  V  X  +  t2,dtldt2  -  j0J0  5lx  +  V 


x  +  t2>dt1dt2 


for  x  >.  0; 


D. 


H0(x  +  tf  x  +  t)  dt  <. 


G(x  +  t,  x  +  t)dt 


for  all  x  >.  0. 

In  all  cases  G(t^,  t2)  is  the  bivariate  exponential  distribution 
proposed  by  Marshall  and  Olkin  (1967)  where 

G(tx,  t2)  =  exp(-  X1t1  -  \2t2  -  X12max(t1,  t2>)  for  t ^  t2  >  0, 
and 


Remark  2.  The  Kiefs jo's  definitions  of  MHNBUE  arc  restricted  to 
the  class  of  life  distributions  for  which  the  following  conditions 

12  : 


hold  among  )i2,  and  )j 
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yl  +  P2  1  V1  +  V2  1 

(i)  — - - - 10,  (ii)  — - -  -—10,  and 


*12 


y12  y2 


(~  -  J-)(Ul2  ' ’*2) 
V“l  V'  "12  ' 


>  0. 


(2.10) 


The  following  example  shows  that  Klefsjo's  definitions  of 
MHNBUE  does  not  imply  our  definitions  of  MHNBUE. 

Example  4 .  Let  X^,  X ^  be  iid  and  X^,  X2  be  HNBUE.  Choose  y^,  y2» 
and  y^2  which  does  not  satisfy  one  of  the  conditions  in  (2.10). 

Then  (X^,  X2)  satisfies  our  definitions,  but  it  does  not  satisfy 
Klefsjo's  definition  of  MHNBUE. 

3.  MHNBUE  closure  properties.  In  this  section,  we  prove  certain 

closure  properties  of  MHNBUE  distribution.  Let  8,,  8~,  B^,  and  8,. 

12  3  4 

denote  the  classes  of  life  distributions  satisfying  the  definitions 
(i) ,  (ii) ,  (iii) ,  and  (iv)  respectively  of  MHNBUE.  Then  we  have 
the  following  theorem. 

Theorem  1 . 

(PI)  If  (Tj, . . .  ,Tm)  eBj  »  any  subset  of  (Tj, . . .  ,Tm)  ^  (1  <.  j  <  4) ; 

(P  2)  If  (T.,...,T  )e8-»  (Tw  . .  .  ,T")  eg  .  and  are  independent,  then 
1  m  3  1  n  3 

(Tlf...,Tm,  T£,...,T')egj  (1  1  j  <  2).  If  (Tx, . . . ,Tr) eg3, 

(T£, . . . ,T') eg^  are  independent  with  identical  distritution,  then 

{Tl'...'Tn'  Tl'-*-'Tn,e83* 

(P  3)  If  (T  , . . . ,T  ) e$  . ,  then  (c.T. , . . . ,c  T  ) eg .  (1  <  j  <  2)  for 
1  m  *  3  11  mm  —  — 

all  c.  >  0  (1  <  i  <  m).  If  (T.,...,T  )g«^,  then  (cT,  , .  .  .  ,cT  )  e.R  • 

1  -  —  1  m  j  1  m  j 

(2  1  j  s.  4)  for  all  c  >  0. 
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Proof. 

(Pi):  This  property  follows  immediately  from  the  definitions. 

(P  2)  :  Let  a  =  (a. , . . . ,a) ,  an  m  vector  and  a'  =  (a' , . . . ,a') ,  an 
'  '  —  1  m  —  in 

n  vector.  Other  vectors  can  be  defined  similarly.  Let 
t 

I  (x)  =  j  (P(T  >  x)/J  P(T  >  yjdy]dx,  where  dx  =  dx^...  fdx^  and 


is  a  multiple  integral  with  m  factors.  I(x')  and  I(x,  x')  are 


defined  similarly.  Then 


m  n 

n  t.  n  t:/i(x,  x') 

i=l  1  j*l  3 


|Ht./I(x)][  t'/Ux')  1 

i=l  1  j=l  3 

—  oc' 

<.  |  P  (T  >  x)  dx  |  P(T^  >  x' )  dx/ 


#o o'  roo 

-  J  j  >  *  *  r  >**" 


) dx'dx 


0'  0 


Similar  proof  works  for  To  Pr°ve  use  the  Jensen  inequality 


(P  3)  :  If  T  e  3^,  c^  >  0  (1  <.  i  <.  m) 


m  r  t  r” 

II  t./  (P  (d  >  x)  /  P  (d  >  y  +  x)dy]dx 

i=1  1  J£  o 

=  II  t./f  [p  (T  >  e)/[  p  (T  >  f)dy]dx 

i=i  1  Jp_  "  -  Jo  -  “  x  " 

2  .» 
i  f  P (T  >  e) dx  =  P(d  >  x) dx . 

Jq  "  ~  ~  J0 
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Here  «  =■  (c^ . emTm)  .  Z  +  *  =  <Yl  +  Xj . Ym  +  *n>  . 

e  *  txi/ci» • *  * ,xn/cm)  and  1  =  (y]/ci* - #ym/e  ) .  Similar  proofs 

work  for  80,  B,,  and  B . . 

2  3  4 

Remark  3.  Since  it  is  known  in  the  univariate  case  that  the  HNBUE 
is  not  closed  under  the  formation  of  coherent  system  (see  Kiefs jo 
(1980) ) ,  the  same  cannot  be  expected  for  the  MHNBUE  under  any  of 
the  four  definitions. 


Remark  4 .  Example  1  shows  that  (P  2)  does  not  hold  for  S^. 

Remark  5.  To  prove  that  MHNBUE  is  closed  under  limits  in  distribu¬ 
tion,  we  need  an  extra  condition  to  guarantee  the  application  of 
the  dominated  convergence  theorem. 


Theorem  2.  Let  {  (T^ , .  .  .  ,T^)  ,  k  _>  1} 

dom  vectors  belonging  to  S .  for  each  k. 

3  st 

weakly  as  k  — **»  and  (T^, . . .  ,T^)  <.  (S^ 
m 


be  a  sequence  of  MHNBUE  ran- 


I£  «•- . Tmk> 


'Ik' 


, . . . ,S  )  for  all  k 
m 


* (T ^, . .  .  , 
■L  kQ  where 


T  ) 
m 


E(  II  S.)  <  *>,  then  (T.,...,T  )  e  B-s  for  each  j. 
i=l  1  1  m  3  ; 


Proof.  Use  the  dominated  convergence  theorem  and  the  appropriate 
definition  of  the  MHNBUE. 

It  is  known  in  the  univariate  case  that  HNBUE  is  closed  under 
convolution.  We  have  not  been  able  to  prove  the  same  result  for 
the  multivariate  case.  Instead,  we  have  proved  the  following  theorem 
which  is  a  special  type  of  convolution. 


Theorem  3,  Let  (a)X  =  (X,,...,X_)  and  Y  =  (Y  Y_)  be  component- 

-  —  I  n  —  i  n 

wise  independent;  (b)  X  and  Y  are  independent;  and  (c)  X  and  Y  c 
for  j  *  1,  2.  Then  W  -  X  +  Y  e  Pj  for  j  =  1,  2. 
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Proof.  Use  (P  2)  and  the  fact  that  X^  +  Y^'s  are  HNBUE. 

Remark  5.  In  Theorem  3  if  X  and  Y  have  the  same  distribution,  then 
W  =  x  +  Y  e  3  . 

4.  Relationship  between  MNBUE  and  MHNBUE. 

In  the  univariate  case,  a  life  distribution  function  is 

said  to  satisfy  the  NBUE  property  if 
•  00  #00 

I  H(x)dx  <  H x ( t )  J  H x ( x )  dx  for  all  t  >.  0. 
t  0 

Rolski  (1975),  and  Klefsjo  (1980)  have  shown  in  the  univariate  case 

that  NBUE  ==*•  HNBUE,  and  that  the  converse  implication  does  not 

hold.  It  is  of  interest  to  know  whether  a  similar  implication  holds 

in  the  multidimensional  case.  Buchanan  and  Singpurwalla  (1977) 

have  given  several  definitions  of  the  MNBUE  based  on  multivariate 

generalization  of  the  univariate  NBUE  distributions. 

We  first  introduce  the  following  definitions  of  MNBUE  as  given 

in  Ghosh  and  Ebrahimi  (1980) .  A  p-variate  distribution  function  H 

;  P 

is  said  to  be  MNBUE  if 


A. 


*  *  *  J  Hp(xy  •  )dxj, . . .  4xp<.  Hp(tj,  •  •  •  ftp) 

‘p‘1 


(Xj, 


,  x  ) dx„  ... , dx 
r  *  P 


for  all  t  ^  0  (1  <  i  1  p)i  and  similar  inequalities  are  assumed 

to  hold  for  all  subsets. 
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B. 

*00  #00  *00  *00 

[•••(  H  ( x  1 9  •  •  •#  X  )  dx  «#»••!  dx  <,  H  (t  # .  .  .  #  t)  f. . .  |  H(x_/***f  x  )  dx  <■/•••  #■  dx 
JJpl  Pi  PP  JjPl  Pi  P 

t  t  0  0 

for  all  t  j>  0,  and  similar  inequalities  are  assumed  to  hold  for  all 
subsets. 

C. 

*00  *00 

H  (x, . . .  ,x)  dx  <.  H  (t , . . .  , t)  H  (x,...,x)  dx 
J  P  P  J  P 

t 

for  all  t  >.  0,  and  similar  inequalities  hold  for  all  subsets. 

Ghosh  and  Ebrahimi  (1980)  have  shown  that  the  following  impli¬ 
cations  hold:  A  =*'  B,  A  C,  and  C  =>  B. 

We  first  explore  the  interrelationship  of  the  A  -  C  definitions 
of  MNBUE  with  MHNBUE-I  to  IV. 

It  is  trivial  to  check  that  MHNBUE-I  =*>  MHNBU-I  =*•  MHNBUE-II, 
MNBUE-II  =*•  MHNBUE-I  1 1 ,  and  MNBUE-III  ;=*  MHNBUE-IV. 

Example  5.  (MNBUE-I  MHNBUE-IV) .  Consider  once  again  example  1. 

Then  (X1#  X2)  is  MMBUE-I.  But  (X^  X2)  is  not  MHNBUE-IV. 

Remark  6 .  The  example  4  shows  that  MNBUE-I I  =£>  MNBUE-IV. 

Remark  7 ♦  Let  X^  and  X2  have  bivariate  Marshal-Olkin  (1967)  exponen¬ 
tial  distribution  with  survival  function 

HjUj^  x2)  =  exp[-  a1x1~^2x2“  *12  max(xi'  x2^ 

x^  >.  0,  x2>.0,  >.  0,  \  >  0,  Aj^2  >  0.  Then  H2(x^,  x2)  in  both 

MHNBUE-IV  and  MHNWUE-IV,  MtINBUE-IJI,  and  MHNWUE-IIT. 
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Remark  8 .  To  see  the  interrelationship  between  MHNBUE  and  other 
classes  of  life  distribution  see  Ghosh  and  Ebrahimi  (1980) .  It  is 
interesting  to  mention  that  multivariate  IFRA  defined  by  Block  and 
Savits  (1980)  implies  MHNBUE-IV. 

The  next  example  shows  that  MHNBUE  class  of  life  distributions 
is  larger  than  the  class  of  MNBUE  life  distributions. 

Example  6 .  Let  X^,  X2  are  independent  and  identically  distributed 
with  the  following  survival  distribution  function, 

.1  0  <.  t  <  1 

_  ?  1  <  t  <  2 

F  (t)  = 

?  2  <  t  <  4 
0  t  >  4 

Then  (X^  X2)  is  MHNBUE- I ,  MHNBUE- I I ,  MHNBUE- I II ,  and  MHNBUE-IV. 

But  (Xx,  X2)  is  not  MNBUE- I  or  MNBUE-II  or  MNBUE-III. 

The  following  example  shows  that  the  MHNWUE  class  is  also  larger 
than  the  MNWUE  class. 


Then  (Xj,  X2)  is  MHNWUE- I ,  MHNWUE-II,  MHNWUE-III,  and  MIINWUE-IV.  But 
(X.,  X2)  is  not  MUWUE-T  or  MNWUE -I I  or  MNWUE-III. 
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